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Abstract 


A semi-empirical  theory  is  presented  for  predicting  the  boundaries, 
velocity  decay  on  the  axis,  and  the  velocity  at  any  point  in  an  air  jet 
or  variable  density  exhausting  into  air  at  rest. 

The  integrated  axial  aicoientuin  equation  is  ^ployed  as  the  basic 
equation  in  the  analysi.s  and  an  approximate  expression  for  density 
variation  is  presented  and  aimlyzed.  The  shear  stress  is  represented 
by  a form  suggested  ly  Roichardt.  This  form  is  siirpler  than  the 
ordinary  Pr^ndtl  expression, 

Calculatioiis  are  made  for  range  of  density  ratios  and  an  analj'-sis 
of  tlia  results  is  attempted, 

A ooTiparison  vd.th  experiment  Is  made  indicatiiig  that  the  ncxirig 
parameter  k is  dependent  upon  the  density  ratio, 

A TT.etho^  for  solution  of  the  mixJ.ng  problem  using  a more  exact 
density  vari^ition  is  indicated  in  Apjriendix 
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THE  MIXING  OF  AN  AHALLY  SYMMETRIC  CCMPRESSIBLE  JET  WITH  QUIESCENT’  AIR 

Bitroduction 

The  problem  of  the  mixing  of  a turbulent  Jet  with  its  i»jLrroundicg 
medium  has  been  analysed  extensively  in  recent  years  because  of  its 
obvious  iaiportance  in  many  engineering  applications  and,  also,  because 
of  its  direct  association  with  the  problem  of  free  turbulence. 

In  1926,  Tollaien(l)  effectively  initiated  work  in  this  field  by 
extending  an  analysis  of  laminar  Jets  to  turbulent  Jets  through  the  use 
of  the  Prandtl  mixing  length  theory(2).  His  theory  was  applied  to  Jets 
issuing  fr«a  a point  oource.  Since  that  time,  many  noteworthy  analytical 
and  experimental,  contributions  tiave  been  made.  Among  these  hs&ve  been 
Kuethr  ■ calculation  (3),  in  193&',  of  the  velocity  profile  in  the  initial 
portion  of  a Jet  exhavisting  into  2 medium  at  rest;  the  suggestion  by 
Reich3rdt(U),  in  19Ul,  of  a new  theory  of  free  tui'bulencej  the  caicula" 
tier!  of  Jot  spread  and  velocity  decay  of  a finite  Jet  codiausking  i.ato 
a moving  stream  hy  Squire  and  Trouncer(3)  in  19UU;  the  experimental 
work  of  Forsuall  and  Shapiro(6),  reported  in  195^0;  and  the  analytical 
and  experimental,  analyses  that  have  been  in  progress  at  the  Engineering 
Expei'iment  Station  of  the  University  of  IlJLinois  since  19iiS»  A large 
peroentvage  of  the  work  that  has  been  dons,  including  those  contributions 
mentioned  above,  deals  only  with  incompressible  flew.  Although  there 
have  been  a ipjja'tjsr  of  experimental  investigations  of  ccss^^reasibla  Jets, 
the  only  tlieoretical  approaches  that  are  t/b11  knc?»n:'i  are  those  of  ?ai(7) 
and  S7:ablew.skl(8), 

Iii  gc-.neral,  one  of  three  approaches  is  u:?od  when  a tbsorstical. 
aruiiysis  of  the  prob.lerr.  is  atte,a.ptedfc  In  the  first  approach,  a soJ.ution 


of  the  turbulent  differential  equations  of  motion  is  sought  in  order  to  determine 
the  velocity  profile  as  a function  of  axial  distance*  Kuethe,  for  the 
incoiipreasible  case,  and  both  Pai  and  Saablewski,  for  the  congsressibla 
case,  followed  this  method,  Kuethe  found  an  approximate  solution  to  his 
equation  using  the  Prandtl  mixing  length  theory  expression  for  the 
txurbulent  shear  sti^ss,  Pai  reduced  the  turbulent  equations  to  a 
solvable  form  after  relating  the  turbulent  fluctuations  to  the  mean 
properties  of  the  flow  through  Taylor's  modified  vorticity  theoiy(2)  and 
representing  the  density  variation  by  a form  of  the  Crocco  energy  integral, 
Ssablewski  assumed  that  a modified  form  of  the  Prandtl  mixing  length 
expression  represented  the  Siiear  stress  and  eliminated  density  from  his 
momentum  equation  by  assuming  that  the  apparent  tui'bulent  viscosity 
coefficient  divided  by  the  density  was  constant  vdfch  respect  to  the 
radial  dijiiension,  A roajor  difficulty  arising  in  tld.s  type  of  treatment 
is  obtaining  a solution  of  the  equations  which  are  always  nonr-iinear, 

A small  perturbation  analysis  ’will  usually  iijieari'-:.e  the  equations  for 
smell  density  and  velocity  vax'iatioiis  and  reduce  the  solution  difficulties. 
However,  large  variations  irr/ariably  necessitate  the  use  of  highly  involved 
numerical  methods  of  solution. 

The  second  approach  is  one  which  uses  an  intf'grated  momentura  equation 
of  the  Karman  type (2),  In  this  type  of  analysis,  velocity  profiles  i>tust 
be  assumed  to  be  similar  throughout  the  Jet  or  at  least  in  portions  that 
are  treated  separately,  and  a particular  velocity  profile  shape  must 
be  assumed  in  order  to  solve  the  res\;lting  equations.  Squire  and 
Troxincer  used  this  method,  representing  the  shear  stress  by  the  Rcandtl 
form,  to  calculate  the  spread  cf  the  .jet  am  the  velocity  decay  along 
®ie  axis  for  the  case  of  incompressible  mixing 


betwp*=!n  a jet  and  a moving  external  stream.  Their  calculations  were 
much  simpler  tha:i  those  arising  from  the  use  of  the  ordinary  differential 
equations,  and  their  assumptions  regarding  the  shape  and  similarity  of 
velocity  profiles  r^ere  rrell  verified  hy  Forstall  and  Shapiro, 

The  third  method  that  has  been  used  is  based  upon  Reichardt*s 
theory  of  free  turbulence.  Essentially,  Reichar-dt  replaces  the  gradient, 
or  conductive,  mechanisia  of  raoznentum  transfer  by  a diffusion  process, 
thus  discarding  the  basis  upon  which  the  Prandtl  expression  for  shear 
stress  is  derived.  Some  experimental  verification  of  Rsichardt’s 
hypothesis  has  been  achievedj  for  example,  Kivnick(9),  at  the  University 
of  Illinois  has  fitted  tlie  data  of  Forstall  and  Shapiro  to  an  analysis 
of  the  Reichardt  type  reaflcmaV-ly  well. 

It  should  be  x'ealissed  that  those  analyses  of  tuirbulent  jet  mixing 
can  only  be  classified  as  seml~enpirical  txeories  because  of  the 
unknc?m  constants  that  always  appear  in  the  results.  These  constaiits 
result  from  the  rather  arbitrary  representation  of  the  turbulent  shear 
stress  and  must  be  evaluated  experimentally. 

It  is  the  object  of  this  paper  to  include  the  effects  of  com" 
pressibiUty  in  the  mixing  problem  by  a method  that  will  allow  solution 


without  a prohibitive  anourit  of  calculation*  The  velocity  of  the  out~ 
side  stream  is  taken  as  zero  to  simplify  the  problem  and,  aI.?o,  because 
the  expe’.'imental  equipment  at  ttds  laboratory  can  be  easily  adapted  to 
such  a condition, 

A method  similar  to  tliat  used  by  Squii  t;  and  Trcioncer  is  employed. 

The  effect  of  density  variation  is  accounted  for  by  ;m  approximation 
that  assirmes  the  density  at  any  point  jui  the  jet  is  a particular  fuiicti.on 
of  the  local  mean  ve3.cclty  and  constant  qn  Anti. ties'  on'iy.  because  it  is 


u 


believed  to  be  not  valid  in  con^ressible  flcfw,  the  usual  Pi-andtl 
representation  of  shear  stress  is  discarded*  A somewhat  simpler 
form,  similar  to  that  suggested  by  Reichardt,  is  used. 

The  transfer  of  mass  and  energy  are  not  treated,  since  experimental 
resxilts,  particularly  those  of  Townsend(lO),  indicate  that  the  mechanism 
of  transport  for  these  quantities  in  a free  turbulence  phenomena  is 
different  from  that  for  momentum.  Therefore,  it  does  not  seem  logical 
to  present  similar  an».lysea  for  the  three  quantities. 
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Theoiy 

/u  Derivation  of  Integral  Equation 

The  integrated  axial  iiwiaenttm  equation(2)  used  in  this  analysis  is 

A brief  derivation  'k113.  aid  in  the  understanding  of  equation  i and  its 
subsequent  use  in  the  solution  of  the  mixing  problem* 

Consider  a circulai-  disk  of  air,  as  sho^m  in  Figure  1,  in  a general 
steady  flow  which  is  symmetric  with  respect  to  the  x axis*  The  change 
of  a.'d.al  momentum  across  the  faces  of  the  disk  per  unit  tine  is 

(2) 

The  change  of  mass  between  the  two  faces  of  the  disk  per  unit  time  is 

and,  since  this  mass  can  only  flow  across  the  r r^s-  siorface,  t)ve 
ari.al  notrientum  change  associated  with  the  mass  flow  change  is 


//I* 


(1) 


ne 


2.?} /I 

force  exerted  on  the  disk  by  the  shear  siu'ess  '2*’  ^ 


(3) 


(U) 

Notice  the  Expression  k ignores  the  change  in  't!  betweeri  z;  and  x g x 
ifhich  is  considered  to  be  a second  order  effect,  Tlv?  pressxtre  g3*a,dients 
in  the  r and  x directions  are  assumed  to  be  negllgibio,  foilCTrlng  the 
usual  boundary  la;rer  and  jet  assurm^tlons,  r.hus  eain'ing  no  uiibalanced 
pressure  force  on  the  disk. 
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A euatnation  of  forces  acting  on  the  disk  of  Figure  1 Expressions 


2j  3,  and  U*  yields 


The  signs  of  the  terms  in  Equation  5 are  dictated  by  the  type  of 


problem  being  Investigated  and  the  form  that  is  to  be  used  to  represent 


the  shear  stress  If,  Fauation  5 can  be  trariSformed  siirply  by 


Lelbnita’s  Rule  to  the  form  given  in  Equation  1, 


Approach  to  Free  Jet  Analysis 


The  derivation  of  solvable  eqiiations  in  terms  of  unknown 


quantities  used  in  this  analysis  is  similar'  to  that  employed  by  Squire 


and  Xrouncer  and.  will  bo  discussod  only  briefly  here.  As  ohotm  in 


Figure  the  jet  is  divided  into  two  regions}  the  core  region,  in 


which  a ’’potential*'  core  of  constant  qiiimtities  exists:  and  the 


developed  regi-on,  in  which  tiie  core  has  disappeared  and  the  properties 


of  the  jut  are  changing  continuous!)'’  •with  r and  x.  Notice  tiiat  the 


mduclng  region  bouixdaries,  x°q  and  are  defined  by  the  spread  of 


velocitys  or  momentusa. 


[t  is  assutaed  then  that  similarity  of  velocity  profiles  exists  in 


both  regions  and  that  these  profiles  can  be  represented  by  p cosine 


vai'iatlon  (see  Figure  3)o. 


Core  region: 


Z ( 


cos  li 


Developed  region; 


(7) 

This  assianption  is  ^vell  verified  by  experiment  in  the  incompressible 
case  (6)  and  early  experiments  at  thiv=?  laboratory (11)  indicate  that  such 
an  assumption  is  still  reasonably  valid  for  sax  extremely  cold  jet* 

Notice  that  both  Equation  6 and  7 become  iaentical  at  tlie  dcrimstreani 
limit  of  the  potential  core. 

New  supipose  that  the  density^  and  the  shear  stress  can  he 
represented  as  functions  of  u,  rc,  and  ri*  Thi.s  will  be  discussed  in 
the  following  sections.  Then.j  in  principle,  the  problem  can  be  solved 
in  the  following  manner; 

Core  Region 

Evaluate  Equation  1 at  r-!:-  Tr. : 


/it 


OL 


Thertifore,  Equation  1 becoraeo 


fa 


sss  o 


or 


,/£«» 


(B) 


Eqoatiori  6 egresses  the  condition  of  conservation  of 


total  Lao!i%’ntiiiuj 


the  ;aoiTiii’ioum  at.  the  jet  exit  vridcii  is 


•will  be  eq^^al  to  the  botal  moiaentujn  integrated  across  the  jet  at  any 
station  x.  Therefore,  the  constant  M in  IJiuation  8 can  be  evaluated 
iiamedistely  as 


» 

Equation  8 can  be  reduced  to  a more  usable  form  by  splitting  the 
integral  into  two  porta 


(9) 


/ * J 

/'t*, 


sr 


(10) 


In  Equation  10 


so  t.hat  Equation  10  become 


•^Aa 


A'€. 


Evaluats  Equation  1 at  rfr  =- 


(11) 


££■* 


therefore.  Equation  1 bscomes 


a 


f^uyi^A  ~ ^ ~ 

at 


(12) 


Ags'ln..  by  splitting  the  integrals,  Fa.uation  12  reduces  to 


.. 


••'  t’- V ■^'  '•  '-1 


ea  ^^^/^  - j^aA^/L 


-hr  "5^  3<v  ^«>2  a 3if 


__ 

According  to  the  assurnption  above  that  and  may  be  expressed  in 
terms  of  n,  ri,  and  Tq  and  that  u is  a knovm  function  of  r.  Equations 
11  and  13  are  recognized  to  be  two  similtaneous  equations  for  the  two 
unknows,  r^(x)  and  rQ(x),  and  in  principle,  can  be  solved  for  these 
unknowns , 

Develoned  Region 

The  derivation  of  equations  in  the  developed  region  is  siitdlar  to 
that  in  the  core  region.  Flcwever,  the  unknovm  Ug(x)  new  replaced  the 
'onki'.cwn  rj^(x). 

Evaluate  Equation  1 at  r-u-  = VqX 


'T''  — O 


y /**  v/2« 


ThpT!  Equation  1 be  comet 


. J"  ^C> 


which  expresses  the  same  relationship  as  Equation  3.  Therefore, 


A ^ 


ar  a.  A, 
51 


Evaluate  Equation  1 at  r-'-  - 


Iv.'i nation  "i  becomes 


— €/i-z 
•£ 


_ ^ ’ '-.f  ^ y/lof-aL 

Sij/’U - |-‘4/ 


(15) 
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Equations  lit  and  1$  are  two  aimioltanfeoua  equations  for  the  two  unknowns, 
rn(x)  and  u^(x). 

It  is  seen  now  that  the  problem  •will  be  completely  solved  when 
ri(x),  and  "'^dx)  are  kiiown,  since  the  velocity  profiles  at  any 
station  x can  be  constructed  by  the  use  of  Equation  6 or  Equation  ?, 

C.  Density  Variation 


As  has  been  iraentioned,  the  use  of  the  two  sets  of  simultaneous 
equations  derived  above  for  the  solution  of  the  problem  depends  upon 
the  possibility  of  expressing  the  density  ^ as  a function  of  not  more 
than  the  three  variables  u,  Tq,  and  The  special  „ase  of  incompressible 
flow,^!=  constant,  has  oesn  created  by  Squire  and  Trouncer.  In  this 
paper,  an  attempt  is  rrade  to  solve  the  mixing  problem  for  the  general 
case  of  a jet  of  uniform  initial  dsnsi-ty  -which  is  ai'bitrarily  different 
from  the  density of  the  outside  rfedlura.  It  is  necessary  then  to 
chocfle  a suitable  f’onrtiori  to  represent  the  variation  cf  density 

Pai  uses  an  expression  of  the  Crocco  type  to  satisfy  the  temperatui'e 
variation  in  the  mixing  region  of  the  free  jet. 


^ su 


Tte  quantities  A and  D are  constants  that  are  deue,nrd.ned  hy  the  bound.2jy 
conditions  of  the  prybleio*  Since  the  pressure  in  the  free  jet  is  assumed 
to  be  constant,  the  density  variation  oo’.ild  be  expressed  as 


(16) 


or,  in  other  words.,  the  density  variation  is  a fianction  of  velocity  alone 


11 
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?.nd  constant  quantities.  It  was  observed  that  the  use  of  Expression  l6 
would  cause  great  difficulty  in  the  evaluation  of  the  integrals  in 
Eqxiations  ?1,  13»  lU?  and  15  when  Fx^uations  6 and  ? were  substituted 
for  the  velocity  variation,  'Kisrefore,  a different  variation  of 
density  is  usedj  however,  its  dependence  upon  velocity  alone  is 
retained.  Appendix  I presents  a possible  method  for  the  inclusion  of 
a variation  of  the  type  siiovni  in  Equation  16, 

The  density  variation  that  was  chosen  is 


(17) 


The  constant  in  Equation  1?  is  determined  from  the  condition  at  the 
nozzle  exit  as 


and,  therefore. 


€1  ~ 

the  density 


at  any  point  of 


the  mixing  region  is  given 


&Z&  u • *-P^ 


(18) 


It  is  realized  that  Squatlon  18  is  orily  ari  approxlTfiatioa, 
of  its  validi^ty  can  be  made  in  tlie  fol3.C5wing  mannere 


Since  M ~ constant,  and 

T fk 


An  estimate 


(19) 
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if  the  special  case  in  wiiich  the  jet  1-ia.s  the  same  stagnation 
enthalpy  as  the  surroxinding  fluid  is  considered,  = T^-j^  , and 


_/V, 


It  is  assumed  that  = Tg  tliroughout  the  nixing  region.  Therefore 
Equation  19  becomes 


« 


or 


r ' a. 


(20) 


Equation  20  is  similar  to  the  ordinaiy  expression  for  the  adiabatic, 
one  dimensional  floir  of  a perfect  gas 


7Si 


/ f' 


r~/  £7l  “ 
a ,!«  * 


except  that  the  constant  Is  substituted  for  the  vari.ablo  (throughout 
the  mixing  region)  a^. 

Consider  a cold  jet;  i,e,,  one  in  which  the  tejiperature  of  the  jet 

is  less  than  the  surrounding  medium.  The  quaiitity  wj.iJ.  always 

be  larger  than  as  one  goes  away  from  the  potential  core  because 

of  the  restriction  yilaced  upon  a^.  The  error  will  be  small  near  the  eos'e 

and  will  Ixicrease  as  the  distance  from  the  core  increases,  Hcwever, 

2 

as  one  goes  away  from  the  core,  u becoires  very  small  and  she  contribution 
of  the  increased  error  in  ' is  believed  to  be  small,. 

for  a out  jet;  is,e.,  one  in  which  the  temperature  of  the  jet  is 
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■c. 


■^3 


Greater  than  the  surrounding  medium;  it  is  impossible  to  have  the  same 
stegnation  enthalpy  in  both  the  jet  and  the  external  fluid,  and  the 
analysis  does  not  apply. 

Suppose  that  ss/j  been  used  instead  of  Equation  17, 

The  equation  that  is  derived  for  the  case  of  constant  total  enthalpy, 
corresponding  to  Equation  20  for  the  variation  is, 


Zs 

r 


/ 


•a-  Uf 


(21) 


u 

In  Equation  21,  the  term‘d  has  a reducing  effect  upon  the  right 

side  of  ih,'  equation.  For  the  cold  jet,  this  is  compensating  compared 

to  the  u'  variation,  no7;ever,  it  should  be  noticed  that  the  compensating 

ef^'ect  for  the  cold  jet  will  pi*ohAbly  be  too  great  at  larger  distarices 

u ^ 

from  the  potential  core  since 

If  a density  variation  of  the 'type 
the  equation  corresponding  to  Equation  20  is 


is  used. 


/ ■» 

~'r 


r ^ •r:;;''  Tr~^  j-t  .e  . 


a"  ^/u. 


(22) 


/ 


term  in  Equat3.on  '22  will  add  to  the  errcr  for  the 


The 
cold  jet« 

Therefore,  it  is  believed  that,  at  least  for  a cold  jet,  the  density 
variation  can  be  represented  by  an  expi’ession  of  the  form 

w 


where  n %,'';u:'les  bstwsen 
ans.i:v5i3:t  the  vales  ol' 


2 anc  3 t'u'ougiiout  the  mixing  region^  For  this 
n 2 was  used.  It  i,&  a.ssixiaed  t>iat,  az  least  for 


I 

! 


111. 


a jet  vrhose  initial  density  is  not  radically  different  from  that  of  the 
surronriding  fluid,  the  density  value  given  by  Equation  8 mil  give  a 
reasonable  representation  of  the  actual,  phenomena. 

It  should  b©  mentioned  that  the  density  change  represented  !:y 
Equation  16  -will  satisfy  the  condition  of  a stream  of  arbitraiy  velocity 
so  long  as  the  assumptions  of  the  analysis  are  not  violated;  for  example, 
the  presence  of  shock  waves  in  the  jet  core  would  spoil  the  assumption  of 
unifonn  velocity. 


Shear  Stress  Bepresantatioo, 


The  representation,  of  the  shear  stress  in  an  arialysis  of  a turbulent 
flow  has  teen,  a controversial  spaqject  fcr  years.  'lbs  sost-  well- 
fci.own  Eethcd  is  that  of  the  renentur  trarisfer  theory  based  upon  the 
Prarvdfcl  Arlag  length  theory  whl.ch^  attempts  to  i«late  the  turbuiesti- 
sbear  stress,  es.&enti*Xiy  & turbulsat  fluctuation  phs-ucasena,  zo  tbs 
ntean  properuxes  ox  tns  fluid  through  a factitious  ”‘.n±xing  length*  (2). 
This  theory  lias  giv-ee  reasonably  good  results  in  imonpressible  pErobl8®iS 
for  deteT5Klning  the  spread  of  ;E33e:3tu.'a,  Kcwirrer,  in.  tb*  C'Orio.adera.ttc3 
cf  fi  t'iirb-ulent  oorrpressifcle  It  -icet  .tot  seem.  Icplca^i  to  use-  scc.h 

a representation,  siiKre  the  cc-src-r^H’jible  turcii.ler?t  equoiions  of  Esott-on 
caiiTiOt  be  derived  cc  give  tas-  sasce  .fera  for  the  turfculeitt.  shear  stress 
as'  the  Incorpresslble  ©cuai-ions* 

F'Cr  this  raas-on,  a sjjepiie  fora  for  the  turbulent  she-ar  stress  i.s 
used  in  thas  paper  instead  of  the  fom. 


^A 


Vi 


1^. 

In  Equation  23>  e can  be  referred  to  as  a txxrbulent  viscosity  coefficient 
and  u is  the  niean  velocity  of  the  flow.  Some  attempt  must  be  made  new 
to  express  e in  terms  of  u,  r^,  and  r^  only  so  th-at  the  equations  derived 
in  Section  ^ ssay  'fcs-  s-olred»  Suppose  tifat  e can  be  broken  into  a dersii^ 
dependent  tens  ajod  & ^ecKKtricaUy  dependent  teiai,  Tihe  staplest  vay  to 
dc  this  Is  tc  say 


^ a 


i ’ 


wr;£TC  £ can  ce  ccasidered  to  be  a turbuien'  ktnesa'r^c  viscoai  y 
coefficient.  It  is  assuted  to  vary  only  (Sithjfas  s^ar^gested  Jr?-  Seichardt  (.b) 
Tfhc  reached  this  cc-ncicifi.ica  through  the  analysis  of  si&ny  incompressible 
to^rissental  data.  Tharefore  ^ is  as^sTjjaed  to  be  proportional  to  the 
T«icith  of  tbs  leixLng  region,  -airhich  is  only  a function  of  x.  The  shear 
stress  Is  then  glvan  by 


(25a) 


D«"/&lopieii  Stagicn; 


"T-FK  /l^ 


C2>b) 


The  -quantitr  1 is  then  experisentaily  cet-ernired  prcportionc'lity 
ccnst-si'jt  of  the  type  i-hat  usualiy  appears  in  ssni-erpiricai  ai:i32:f'ses 
3T^  has  the  diirenssicns  of  a velocity. 
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£»  Soltttios  or  the  Eq’iiatxcgts 


Core  aagi' 


It  is  nan  possible  tc  express  Efcj'iiatiosis  11  axd  15  iri.  teras  oi  oil;/ 
tbe  t?«>  depezsient  Tariables  srid  Sabstitsticr.  cf  Iqaatiofl  iS 


Tor ^ sad.  Squaticn  25  a fsr^giires 


a.  V J 


(25  i 


^ »ncX 


>w<vg* 


^ I J 

be  rirht  sice  cf  equation  2?  is  obtained  izr  evaluating 


i 

J 


r..  + 

Sltati4,>.  25a  st  r = — :—_ — i Tnisa  Eqpmtian  6, 


Ilw  Inteiprals  in  FniiatiOTis  26  and  2?  can  be  eraljuated  by  ordinary 
jas'thods  after  'so'risti'fcatiori  fcr  j.  'rrf  .E^juation  6.  Tee  evsi^sxed  .inbegr;als 
tre  both  equations  can  be  seer,  to  Involve  cnlv  the  variable's  r,^  ai>i  r^, 
Trcn  Tedocti-COj  the  ncryKiiri3n.slo;<jal  ecuations  ■i*E*'ror:e 


^ --f {£-*)- p-J 


{ Oir: 


(29) 


,?■'■  r--¥?  ' ; 


17, 


o i 'jC 

where  Rr,  = — , Ri=’3»r>  X = 77*  The  lettered  constants  are  fid.\-en  in 
ri  1 * 1 

Apcendix  II.  Bjuations  28  and  29  are  similar  in  form  to  those  derived 
by  Sqtiire  scrd  Trouncer  vfho  give  an  intricate  method  of  solution  for 
Rq  and  R^.  This  method  is  not  follovred,  however,  since  it  was  felt  that 
a graphical  integration  would  yield  results  more  quickly. 

The  method  used  to  solve  Equations  28  auid  29  is  as  follows: 

Equation  28  was  solved  for  Rq  in  terms  of  Rj_  and  also  differentiated 

, , 4.  jL  -y-  4.^  ■ I • .t  « 

ivlth  respect  to  2 td  gave  . Rq,  , and were  suostiuuteo 

into  Equation  29  and,  upon  reduction,  an  expression  relating  and  X 
was  foixnd, 


(30) 


where  t is  a function  cf  and  {See  Appendix  II), 

Using  the  boundary  conditica  at  the  Jet  exit 

% = 1 at  X = O 

Iquati'csri  jO  carj  be  inte:grate£i  to  five 


(31) 


The  fUiiotion  f was'  c-alOTlated  for  values  cf  var-j-ing  between  1 

at  tb*  rcszsle  scd.t  arsd  O at  the  eid  cf  the  pctential  cere  aj;<d  for  various 

'deaslty  ratio  values,  2, 1^25  ^ .166'?,  yalaes  ef  X 

.cs.srrespersiir.g  to  toe  chosen  E..  and  vslaes  wisre  then  calculated 

sZ’S’G'ir^iLc  ally  arai'  R sfas  calc-letad  fna  Equ-ation  28. 


'i-'ure  h ■..•voffs  iha  calcrbLatei  and  v.iluss  ss  functicus  c-f 

o _ 


ran-ge. 


.eaves  Tne  _erigT-'.  of  th.s 


P/ 

a i'iu-xitlcn,  cf  de;:-s:.ip'  ratio 


IS, 


Be?gl0p8d  BggLcgut 

ScbstltotioKi  a?  $qwa.tic«a  18  far  ^ssA  B^aatloa  2iib  fcT'  '2T  into 
Equatioeis  Hi  and  15  giTOS 


(32) 


snd 


5^ 


-/~^T  ^J. 


^ C£<- 


Trae  ir..te,gr*l”  'BTaa:,icit\s  32  arsd  33  again  can  be  eTSiiiated  urb'ea 
sabstl  tctios  Xcr  is  mde  nith.  Sqiianiaa  7*  I3»«  n«i!-iiiJ5ertai'Ciial  eqiiatioiis 


3*:CSlf8 


(3li> 


(33) 


/V/^  \ ^ ^ J’C  -^r 


wnere  ~~  j C — 


'Tlie  Iet-b3re4  coTi?tant^J  are  given  In 


^ * y€* 

Appendix  II,  SabstlPuticn  into  Bjaaticn  far  frcu  Equation  3*J  .rss;jiits 


in  aa  ea^irtts-aioa  involviag  onlj'  IF  and  Xj 

^ ^ -r 


(36 


jS? 

rtiers  g (??'  is  s.  fuRvct-xm  ctf  5 and  fmiy  (See  Appendi 
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boundaiy  condition  that 

U = 1 at  X --  X» 

nrhere  is  ths  length,  of  the  pot-eirtial  ccre>  is  then  enplcyed  to 
integrate  Sqaatlon  3^' 


The  .fijnction  g (U'  ■ras-  cai-cila^ec  for  rslues  of  I'  and  in  tixg' 
following  ranges: 

l.  X > U > .12$ 

2.1*2$  > > 1.667 

Valt!«3  of  — X correstJoautB;?  to  U caleolated  ^atihically  according 

€JLj> 

to  equation  37.  values  larere  Cclc::ian«d  fro©  t-^aatioR  3h. 

The  '.'ariatio.i  of  axiai  velocity  is  piot'ted  in  ligare  6 for 
cifiersnt  'valses  of  . Flgua-e  ? sr»Tr£  tbs  spread  of  the  jet  as 

^ 1 y .*-  ^ j.  A Z ... 

— A — *■  A -ie  .***%-^  * 


07) 


ii; 


2D, 


» 

5 

i 

I 


Solution  of  a IFto^viiycilar  Mixing  Prcfolea; 

?lgQ793  8 and  9 present  crosspljsts  of  Figures  U,  6^  and  7 shewing 
the  spread  of  the  Jet  .sud  the  decay  of  velocity  along  the  axis  respectively 
as  a function  of  4}2  fer  various  X values.  %hen  it  is  desired  tc 
pre:iict  tfcs  cha’-acteristics  cX  a Jet  cf  knewr.  oenslty,  or  teaperatuire, 
assusing  the  presstare  to  be  coisttanij  exhaustlag  into  a still  Taedltaa  of 
ioiowc  aeasitw.  Figures  5»  5,  and  9 are  easily  used  to  find  the  pt)t«tti,al 
core  length,  ths  Jet  bcuadarLea,  and  the  velocity  decay  on  the  axis* 
SquatlOK  6 or  'Squotices  7,  the  cosine  Telocity  variations,  -51111  then  give 
the  Telocity  at  any  point  In  the  J«t» 

It  is  the  purpose  of  this  paper  to  discas-s  the  effect  of  densi-|gr 
'jT/m  a Jet  discharging  into  s-tilX  air  -only,  Hcsserrir,  it  seeas  that 
the  case  cf  an  external  strea®  -veloci-by  not  eoual  to  sero  could  be 
taken  into  acenmt  quits  easiiv  by  Ihe  cethod  gtrea  bjr  Kucheai&£in  and 
Weber  (12!  -a  chan  re  of  th.e  ccerdiiiatE  systen,. 

hffect  of  lisnsity  '/ariatlon  Upcai  Jet  Charac  tori  sties; 


Q-jaiitatlvely,  the  effact  of  vsrylng  density  upon  che  characteid-Stics 
esf  a jet  spresdiTiig,  into  a free  str-eac,  nX  aero  veiocirj  can  be  seen  in 
Fi’garea  6 and  ?,  It  is  nctioei  that  the  pcteirtial  core  becciaes  IcTiger 
as  the  density  ratio  increases  (Figure  5)*  A hot-  jet,  lew  Taiuesj. 

not  CEiIy  iitas  a s'r’C-rter  pcteat-ial  core,  but  its  centerline  velocity  also 
■drops  eff  lEUcii  faster  -sri-feh  axial  distance  than  a cell  jet*  I'his  is  lu 
geTTfferal  agre'ea.ein.i  ici'-th  C.crrsir,  isiii  iTberjitij'  «rr  otsera-ed  that  a '•isarE. 
jet  ST-rsads  faster  and.  y-harefcicp.  ite  certarLuie  -veli-'c-ity  dfe-cr-i-7as€-3  isrore 


- 
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rar-idly  th&c  a cool  jet, 

Ths  spread  of  tJie  jet  boandarles^  and  is  iirterestlng  vdbexi 

0f 

it  is  noticed  that  ratil.  an  .asdal  distance  of  about  X = ,05l5  is 

rs'acbed  the  hottest  jet,  «l66?,  has  the  least  outward  sjjread, 

dr  . 

Its,  iasrard  spread  is  the  wsst  rapid,  Hcnrever,  after  ^ X = ,3o7 
this  Jet  has  the  rrost  rarsid  tratmsrd  spread.  These  effects  can  be 
esxjrlained,  qjialitativeljj  as  follows: 

Consider-  a hot  and  a cold  jet  of  equal  Initial  Telocity  eshaustizsg 
frras  noasles  of  equal  radii  into  sarrotrndEijigs  of  equal  density.  The 
total  s!C3®entaffl  flujc  of  a jet,  nmch  will  te  pj-oportional  t?111 

*06'  less  for  tloe  hot  thar!  for  toe  cold  jet.  For  tiais,  reaace.,  a c-oi3s>aritiTe'l,7 
srall  asjo'unt  of  Boneatun  •cransf.ar  Tfill  be  needed  in  the  hot  Jet  'case  tO' 
destroy  the  potential  cere  cor^yletsly.  Cansequently,  a ssatller  arsrmt- 
of  outside  sir  is  needed  to  be  entrained,  causing  less  outsard  spread 
fc?r  the  het  Jet  in  the  core  region*  iUTter  ti»  core  has  disappeared, 
t>*»  vpjiw’ty  on  the  axis  decays  wore  rapidly  in  the  hot  jet,  also  because 
of  the  reasc/i  just  given.  Hiwever,  the  ccHisenration.  of  ECE»entiHE 
conditicri  sscst  'be  pfreser5?s.<i  tta'is  causiag  the  crater  boaiydary  of  the  hot 
jet  to  spread  rapidly  after  the  cere  is  destroyed, 

'TJve  MLxinr  Farasjster,  _h_- 

U-; 


fJritil  now,  c'te  results  oi  the  c-aicuiatioas  have  'begn  present'ad  .in 
ter?^  d ~ ^ & teieg  the  propcrtlonallty  constant  in  tls  shear  stres,s 

s-epresenca-tioti,  Eq-aaticn  2^.  first,  it  appears  t&at  the  absolute 


velccit-y  of  t:5e  jet  cxhaiEstiaa  frcsi 
infloe’ice  uptKi  the  axial  oiciensicn. 


the  nossle,  b ^ , has  an  iTcpcrt^rft 
neyever,  esperiseriai  results  f 12'* 


show  quite  clearly  that  the  characteristics  of  the  Jet  depend  primarily’- 
upon  the  ratio  of  the  extrenal  stream  velocity  to  the  Jet  exit  velocity 
and  only  negligibly,  if  at  all,  upon  the  absolute  velocity  of  either. 
Therefore,  since  this  analysis  was  derived  for  the  free  jet  casej  i.e., 
for  the  particular  velocity  ratio  of  aero,  neither  the  velocity  nor  any 
velocity  ratio  should  enter  into  the  determination  of  the  jet  character- 

jf 

istics.  For  this  reason,  it  is  believed  that  the  quantity  should 
be  considered  as  the  eiapirical  constant  instead  of  K alone, 

A value  of  was  doterained  frois  the  data  of  Pab3t(lii),  which 
was  taken  at  a density  ratio  ~ A value  of  ~ = .0038 

was  found  to  match  the  length  of  the  potential  core,  which  he  found  to 
be  12*0  nozzle  radii*  A comparison  between  Pabst’s  velocity  decay  data 
and  the  velocity  decay  predicted  by  this  ana.Lv3is  (Figure  9)  using 
—T’  = ,0038  is  given  in  Figure  10*  The  agreement  is  reasonably  good 

up  to  the  limit  of  the  experiment  TrMch  is  about  25  nozzle  diameters 


dow»ivStresm 

/C 

There  is,  however,  an  indication  tiiat  varies  with  density 

ratio  ^/jS-  Data  talcen  by  Corrsin(l5)  at  a density  ratio = ,965 
result  in  a value  of  *0059  when  the  theory  is  matched  vdth  experiment 

sat  the  end  of  the  potential  core* 


The  re.sult.s  of  early  experiments  at  this  laboratory  (11)  aiso  shovvn 
Figure  11  for  the  case  of  an  extremely  cold  jet,  = 2,30,  When 

evaluated  at  the  end  of  the  potentxal  core,  as  defined  by  a .sudden  change 


in  center  line  velocity. 


- »C033»  The  theoretical  agreement  '.dth 


velocity  de.,.Ay  downstream  is  then  fairly  good;  no’wever,  the  number  of 
experimenual  data  are  riot  siiff;.cier..t  for  an  accur-at'  eval/eation.  It 
sheuia  be  noted  that  the  jet  exit  Mach  numbo:*  was  aporoxirin.’.eiY  2.6  in 


these  exp'eriinsMts  and  weak  shock  disturbances  present  i.n  the  irdti.al 
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portion  of  the  Jet  mi"ht  have  changed  the  length  of  the  potential  core 

from  its  ideal  value.  It  is  believed,  though,  that  shock  effects  were 

Jt 

s;rall  and  that  the  value  of  = ,0033  is  reasonably  accurate. 

Therefore,  there  appears  to  be  a variation  in  the  mixing  parameter 
ti  \7ith  density  ratio,  which  was  not  accounted  for  in  the  development 

£i4 

of  the  theory  being  presented,  A systematic  experimental  deterndnation 
of  will  be  a part  of  future  tests  conducted  at  tbas  laboratory. 
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Apperu^lx  I Solutdon  of  the  ItLicing  Problem  Using  a Linear  Velocity 
Profile  and  the  Crocco  Integral  for  Density  Variation 

Equation  l6  is  an  expression  of  the  Crocco  integral  type  and  is 
an  exact  solution  of  the  general  energy  equation  in  the  Pi’andtl  nuniber 
or  the  flcnv  is  assumed  to  be  iinity  and  if  the  pressure  is  assumed  to  be 
constant.  It  •wcruld  be  desirable  to  represent  the  density  variation  in 
a jet  by  such  an  expression  instead  of  the  approxiination  given  by 
Equation  18,  Equation  16  can  be  rewritten  as 
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y/here  a<  and  can  be  deteimned  by  the  boundary  conditions  of  the 
given  Jet. 

The  applicability  of  the  method  of  solution  presented  in  this  paper 
depends  upon  the  possibl3.ity  of  being  able  to  express  integrals  of  the 
form 


(I 


ao  algebraic  functions  of  their  limits,  since  the  limits  always  contain 
at  least  one  of  the  unknowns  of  the  problem.  Tliat  is,  expressions 
similai'  to  Equation  1.2  must  be  Integrable,  The  assumption  of  the  cosine 
shaped  velocity  profile  eliminates  the  use  Equation  1,1  as  a density 
vai'iation  because  of  the  resulting  complicated  fonn  of  the  integrand, 
in  Equation  1,2  when  substitution  for  u is  made, 

;m  iiitcre,slving  calculation  lias  teen  made  by  Fitlcin(ll)  at  t.iiis 

He  etp'loyed  a lirieai’  velocity  variation  and  ’c-eealculated 


J.abci'atory, 


the  DJicorapressible  problem  of  Squire  and  Troimcer  for  the  one  case  of 
no  external  velocity,  obtaining  results  almost  identical  vri.th  theirs. 

This  shows  in  pajrt  that  the  aasiuned  velocity  profile  is  of  little 
importance  in  the  method  of  solution,  a result  typical  of  analyses 
employing  integral  equations.  The  velocity  profiles  used  by  Piticin 
are 

Core  Regions 

(1.3) 


(I.ii) 


^ * C /I*  *"/?*/ 


Developed  Regions 


ic-a.  ^ 


Now,  if  it  is  assumed  that  linear  profiles  can  be  used  in  the 
compressible  let  case  without  causing  an  appreciable  change  in  results 
(a  method  of  checking  this  assumption  would  be  to  repeat  the  analy^sis 
presented  in  this  paper  using  linear  profiles  and  compare  calculated 
resxilts),  the  problem  may  be  solved  using  the  more  exact  density 
expression.  Equation  T.l.  This  can  be  seen  by  substituting  Equations 
1,1  and  1.3  into  Equation  1,2  which  gives 


(1.5) 


Upon  reduction,  Eauation  1,5  Wisy  be  split  into  3 integrals,  the  most 


co:npLj.cated  being  oi  the  form 
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Equation  1.6  is  rational  and  may  be  evaluated  easily although  tediously. 


Append! K II 


Constants  and  Vi-inctions  Appearing  in  Equations 


Equation  2S): 
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Equation  29): 
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Equation  30) 
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Note;  loTser  case  letters  refer  to  constaiits  listed  for  EX:;ua.tiDns  28  and  29 
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Equation  3ll): 
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Equation  35); 
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NOTiCS;  GOVERNMENT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  OR  OTHER  DATA  1 

ARE  USED  FOR  s\NY  PURPOSE  OTHER  THAN  IN  CONNECTION  WTHI  A DEFINITEI.Y  RELATEfj 
COVEPN^IEHT  PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMEN’r  THEREBY  INCURS  I 
NO  RESPONST.BXl.iTY,  NOR  AIR/  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERMMSNl  felAY  RAVE  ■:‘'OR?yrULATSD,  FURNISHED,  OR  IN  ANY  WAY  SUPPI.IEO  TRF- 
SAID  DRAPTNGS,  SPECIFICATIONS,  OR  OIHER  DATA  IS  NOT  TO  BE  REGAiU^SD  BY 
IA!.1Y1ICATICK  OS,  OTHERWISE  AS  IN  AK£  MANNER  LICENSmG  THE  HOLDER  OR  ANY  OTRER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PER^ilSSION  TO  MANUFACTUF^ 
USE  OR  SELL  A?r/  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  itELATED  THERETOj 
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